Electron State Entanglement as Double Ionization Result 
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A theoretical model allowing to explain the behavior of a correlation function of electron pair in 
process of double-ionization has been built. It turnes out the correlation functions' behavior can 
be explained by non-orthogonality of entangled wave functions of electrons and the change of spin 
state of a part of electron pairs in the process of atom double ionization from singlet to triplet. It 
is shown that properties of correlation function are in the best accordance with experimental data 
under assumption that the probability of electrons' transition into triplet state is of the order of 
magnitude 0.5 ore higher. 
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PACS numbers: 03.65.Ud, 34.80.Pa 



I. INTRODUCTION 

Effect of double ionization has been described recently 
jl], [| § . A characteristic feature of this phenomenon 
is the simultaneous release of two electrons accompanied 
by correlation of their momenta. The dependence of the 
correlation function on the difference between the mo- 
menta of the electrons is given in for several series of 
experiments. 

A significant error in measurements of the correlation 
function dependence on the difference between registered 
momenta allows to fit some theoretical models to the ex- 
perimental data obtained. Some models of double ion- 
ization partially explaining the correlation of momenta 
of the released electrons were considered in papers [1,5], 
however the initial and final states of the electron pair 
were assumed to be singlet. In contrast to those reports 
we assume that some electron pairs undergo a singlet- 
triplet transition in the double ionization process. 

A substantial correlation of the electron momenta 
found in the experiments indicates essential modification 
of the electron states, since, at least in helium, the elec- 
trons are in a singlet state and have identical momentum 
wave functions and therefore the calculated correlation of 
momenta is quite weak [1]. The entanglement of the mo- 
mentum parts of electron wave functions in double ioniza- 
tion causing essential correlation of the monenta, can be 
explained, in our point of view, by a singlet-triplet tran- 
sition in a part of electron pairs. Magnetic field created 
by the quickly moving charged particles (first of all, by 
nuclei), can change the spin state. Electrons in discrete 
states in the Coulomb field of a nucleus cannot make 
a transition to a triplet state, since it would lead to a 
substantial increase in energy (ls-2p transition), whereas 
there are no such restrictions for transitions between con- 
tinuous states. 

In this report we have shown that the non- 
orthogonality of entangled electron wave functions and 



the change of a spin state of electrons substantially in- 
fluence the properties of the correlation function. By ne- 
glecting the Coulomb repulsion of released electrons and, 
as a consequence, change of the correlation of momenta 
after ionization, we suppose that the electrons form an 
entangled pair either in spin (singlet state) or momentum 
(triplet state) space and compare the calculated correla- 
tion function with the results of experiments. Theory 
and experiment agree when singlet and triplet states are 
equally probable. 



II. CALCULATION OF CORRELATION 
FUNCTION. 

Choose the wave functions of released electrons as 
<Pi,2 (pi,2,P,P,*) = 



(27R7 



exp 



V 



4ct 2 



where P is the total average momentum of electrons, p 
their average momentum relative to the center of mass 
and <7 is an uncertainty in momentum of each electron. 
Average momenta of electrons in these states are 



(Pl,2) 



P±p 



(2) 



Average coordinates of particles correspond to free mo- 
tion 



ri,2 



(P±P)* 

2m 



(3) 



An uncertainty in coordinate of each electron at the ini- 
tial moment is equal to ^ and increases if the electrons 
are assumed to be freely moving 
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At the initial moment the electrons cannot be found 
at too large distances from the location of atom before 
ionization and, therefore, the uncertainty of coordinate 
is restricted. Thus, the uncertainty in momentum is re- 
stricted from below due to spatial localization of initial 



electron states. 

Taking into account the change of spin state of elec- 
trons, the spatial part of two-particle wave function in 
momentum representation can be either symmetric or 
anti-symmetric 



*S=o (pi,p 2 ;P,P,t) = 



<pi (pi; P,p, t) ¥2 (£2; P,p, t) + <p 2 (pi; P,p, t) <pi (£2; P,p, t) 



(5) 



*s=i (pi,P2;P,p,i) 



<£l (Pli P,P, t) <f 2 (P2i P,P, t) - i£2 (Pli P,P, t) <£i (g2j P,P, t) 
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These functions are normalized and satisfy the 
Schrodinger non-stationary equation 



9*5 

iH i>r = H * 8 > 



(7) 



where Hamiltonian describes a pair of free particles H 
P1+P2 

2m 



An overlap integral for tpiand ip2 is 

J = J <p* (p; p P, t) <f2 (p; p,p, t) dp = 

= ex P (-^) (8) 

and depends on average momentum of particles relative 
to the center of mass only. 

Let us consider the statistical characteristics of double 
ionization process. We assume that after ionization the 
electrons gain a total average momentum P , the differ- 
ence between average momenta being p and total spin 
of system S . We denote the differential cross-section of 
this process as 



e s (p,p) = 



das (P,p) 
dPdp 



(9) 



Average characteristics of the state are P , p and S 
. Their estimate can be made taking into account the 
statistical results of measurements. 

Then the electrons freely move towards detectors and 
are registered with some momenta p! and p 2 that are 
not equal to average electron momenta and can be 
very different from the average values. For each specific 
state of an electron pair the differential ionization cross- 
section $g (pi) P2i Pi Pi t) with registration of fixed mo- 
menta pi and P2 through the channel with average elec- 
trons' momenta (P1.2) = and total spin S is equal to 
the product of the differential cross-section of ionization 
into this state and probability of registering momenta pi 



and P2 in this state. This probability is determined by 
the square of two-particle wave function modulus, so the 
differential ionization cross-section through the channel 
with average momenta ^2 , total spin S and registered 
momenta pi and p 2 is 



$s (pi,p 2 ;P,p,t) = 
e 5 (P,p)|*s(Pi,p 2 ;P,P,t)| 2 . 



(10) 



A differential ionization cross-section with registering 
fixed momenta pi and p 2 is determined by the sum over 
all possible states and equals to the sum of integrals of 
products of 6s (P, p) and the squared modulus of a two- 
particle wave function \^fs (pi, P2; P, p, t)\ 2 , this sum is 
taken over spin states: 



$(Pl,P2) = J! 



das (Pi,P 2 ) 
dpidp 2 



= 2 / e s (P, P )|* s ( Pl ,p2;P,p,t)| 2 dPdp. (11) 
s J 

Now we can determine the cross-section of ionization 
with registration of a momentum p for one of the two 
electrons and the total cross-section for double ionization: 

p(p)= J $(p,p')dp', (12) 

Ntot = j $(p,p')dpdp' = 

= J2[ e 5 (P,p)|* s (p,p';P,p)| 2 dP#dpdp'- 

s J 

= f (P, P) dPdp = N + N U (13) 
s J 

where N = J 6 (P, p) dPdp and N x = 
J 0i (P, p)tiPdp are the total ionization cross-sections 
into singlet and triplet states respectively This gives 
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us a possibility of obtaining measurable characteristics 
intensity of registration of two electrons emitted in the 
same act of double ionization 

I cor (Ap) = (Ap) 2 J $ (pi, pi + Ap) d Pl dn (14) 

along with intensity of registration of two electrons ir- 
respective of the fact whether they were emitted in the 
same or different acts of double ionization 



(Ap) 



(Apy 

N tnt 



/p(Pi)p(Pi + Ap)rfpidn, (15) 



that determine the correlation function 



R(Ap) 



(Ap) 



(Ap) 



1, 



(16) 



where Ap = |pi — P2I is the difference of momenta of 
registered electrons. It is this function that is observed 
in experiments on double ionization. 

The differential ionization cross-section $ (pi,P2) and 
the intensities I cor and I unC or depend not only on two- 
particle wave function, but on the probability of ion- 
ization through a concrete state of the electron pair as 
well, and this fact is represented by the factor 0$ (P, p) 
in (11 13,n4,n9). If the only bossible state of electrons 



with certain average momenta and fixed spin is formed 
in the process of ionization then the integration over 
states disappears. Therefore this factor is replaced by 
the number of registered events of double ionization N to t 
and the correlation function is determined by a two- 
particle wave function in this state. If a pair of elec- 
trons with the same probabilities Ps (P,p) undergoes a 
transition to a state with total spin S and average mo- 
menta (pi.,2) = P ^ in the process of ionization, the 
differential ionization cross-section though each of those 
states is equal to 65 (P,p) = N to tPs (P>P)- As a con- 
sequence, the correlation function does not depend only 
on the properties of a wave function of single state of an 
electron pair - it is a characteristic of ionization process 
average over all possible states. 

The most essential experimental data can be summa- 
rized as follows. Firstly, alternation of sings (negative 
for a small difference of momenta, positive for momenta 
comparable with the atomic unit, and negative for large 
differences of momenta); secondly, a maximum lies be- 
tween 2 and 3 a.u. 

In the general expression for differential ionization 
cross-section obtained above there are no assumptions 
about character of electron motion from a point of ion- 
ization to detectors. If we assume that the electrons 
move freely, the following simplifications can be intro- 
duced. First of all, in the momentum representation 
the square of modulus of a two-particle wave function is 
time-independent and the correlation function does not 
depend on the time interval necessary for the electrons 
to cover the distance from a point of ionization to de- 
tectors. With the aim of further simplifications in the 



expression (|llj) let's consider the fact that according to 
(|l|) the squared modulus of a two-particle wave function 
is almost independent on average momenta P and p as 
long as their values are smaller then momentum uncer- 
tainty of each electron a . By assuming that after ion- 
ization the electrons have such an average momenta we 
can replace the differential ionization cross-section with 
registering electrons' momenta pi and P2 by expression 



(17) 



where integration over momenta of resulting states gives 
only the total ionization cross-sections with spins S = 
and S = 1 . Let's denote the probability of a singlet- 
triplet transition as / . In accordance with (|l^) the 
differential ionization cross-sections can be written as 
N = (l-f)N tot , Ni = fN tot . A differential ion- 
ization cross-section with registering of fixed momenta of 
electrons 



N,, 



$(Pl,P2) = 

(/ 1*1 (Pi,P2,t)| 2 + (1 - /) |*o (Pi,P2,i)| 2 ) (18) 



is a product of the total cross-section and sum of prob- 
abilities of registering electrons with fixed momenta in 
singlet and triplet states. Thus, if we assume that aver- 
age momenta of electrons are small in comparison with 
uncertainties of their momenta, the probability of ion- 
ization through a triplet state will be determined by the 
product of a triplet wave function and probability of spin 
flip, /. Then the uncertainty of momentum of electrons 
a acts as a scaling factor, and the shape of correlation 
function is mainly determined by the ratio of probabili- 
ties of ionization producing entangled pairs in singlet and 
triplet states. 
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FIG. 1: Correlation functions for processes of ionization into 
states with different probabilities of transition to triplet state. 

A correlation function for a process of ionization into 
states with the same values of momentum uncertainty 
and different probabilities of transition to a triplet state 
is shown in Fig. [I]. These graphs show that the cor- 
relation of momenta decreases along with the probabil- 
ity of transition into a triplet state. It is clear that 
the spatial part of a two-particle wave function will be 
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a product of one-particle wave functions if average mo- 
menta of electrons are identical. As a consequence, the 
squared modulus of a two-particle wave function is equal 
to |*o (Pi,P2,£)| 2 = \ip{pi,t)\ 2 |^(p 2 ,i)| 2 corresponding 
to an independent distribution of two random values. 

Now we consider the dependence of the correlation 
function on the average momenta of electrons. We as- 
sume that the value of 0s differs sufficiently from zero 
in some area near average momenta P, p, and the width 
of this area is noticeably smaller than the uncertainty 
in momentum a. In such an area a two-particle wave 
function ^5 is essentially constant and the differential 
cross-section of ionization with registering of momenta 



Pi and P2 of electrons can approximately be replaced by 

$(pi,p 2 ) = ^7V s |*s(pi,p 2 ;P,p,<)| 2 , (19) 

s 

where N$ — J 6s (P, p) dPdp again stands for the ion- 
ization cross-sections through the channel with spin S 

Under these assumptions the properties of the corre- 
lation function R (Ap, p) depend on the difference of av- 
erage momenta of electrons p and do not depend on the 
total average momentum P. 

The intensities will depend on p as 



exp I — I J' I sinh z 1 2 sinh z ^ 



>** - \ ' K 1 -/) ^T-n 1 + I (20) 



2y^a 3 1' 1 + J 2 I - J 2 



N tot (Ap) 2 \ II ' (A + B + C), (21) 



cxp(-^)j 2 



A = (i-fY 



B = f 



1 + 2 J 4 + J 



2 sinh z _|_ g j5/2sinh_z/2 



1 + 2J 4 + j2Sinh£ _ 8J5/22SW2 



C = 2/(l-/) 



1 - 2 j 4 



j 



2 sinh z 



(1 + j 2 ) (1 - j 2 ) 



Now we consider the dependence of correlation func- 
tion on p separately for a singet and triplet state. The 
correlation function for each of the states depends on ex- 
pression z = 



and the overlap integral dq) 
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FIG. 2: Correlation functions for several values of difference 
of momenta for electrons in triplet state. 



Ro(A P ,p) 



2 (aSsp + !)(! + J 2 ) 



1 + 2 J 4 + J 2 



sinh z 



gj 5 / 2 sinhz/2 



2( 



sinh ; 



1) (1 - J 2 ) 



1 + 2J 4 + J 2 sinhz - 8J5/2 sinh2 / 2 



(22) 



(23) 



Graphs of correlation function for several values of p 
in a triplet state are given in Fig. |[ They show that 
an increase in p affects only the character of correlation 
but the state itself remains strongly correlated. Starting 
from p ~ a the correlation function loses its second null 
and remains positive for large Ap, which contradicts the 
experimental data [1]. 



Dependence of correlation function on p for a singlet 
state is shown in Fig. [| In this state, on the contrary, 
the correlation increases along with the difference of aver- 
age momenta as the influence of symmetrization of one- 
particle wave functions becomes noticeable. Neverthe- 
less, the correlation function has only a single null. 

These graphs indicate that properties of correlation 
function of singlet and triplet states are different. Such 
distinctions are most visible for large differences in the 
momenta registered. When the difference of average mo- 
nenta of electrons is small (p <C <r) , the correlation func- 
tion is negative only for a triplet state. In the area of 
small differences of momenta registered Ap ~ la.M. the 
correlation functions of a singlet and triplet state are neg- 
ative. For a triplet sate the correlation function is equal 
to -1 at zero, whereas in singlet state it depends on p and 
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FIG. 3: Correlation functions for several values of difference 
of momenta for electrons in singlet state. 



can tend to zero. 



III. COMPARISON TO EXPERIMENT. 
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FIG. 5: Approximation of experimental data for ionization of 
He by C 6+ 

The uncertainty of electrons' momenta produced af- 
ter double ionization should depend on the process of 
ionization. This uncertainty can be determined from ex- 
perimental results. 

The least square approximation of the experimental 
data given in [1] leads to an uncertainty of 0.22 a.u. and 
error of approximation of 15% for ionization of He by 
Au 53+ (Fig. ||), uncertainty of 0.55 a.u. and error of ap- 
proxi mation of 19% for ionization of He by C 6+ (Fi g. ||), 
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FIG. 6: Approximation of experimental data for ionization of 
Ne byAu 53+ . 



and uncertainty of 0.39 a.u. and error of approximation 
of 18% for ionization of Ne by Au+ 53 (Fig. §). 



FIG. 4: Approximation of experimental data for ionization of 
He by Au 53+ . 
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IV. CONCLUSIONS. 

We have considered the influence of a spin state change 
of an electron pair in process of double ionization on 
the correlation function properties. Taking into account 
a large scatter of experimental data-points, the results 
of our approximation show that experiments on double 
ionization are in a good agreement with the assumption 
about the change of the spin state of a part of the elec- 
tron pairs. The negative sign of the correlation function 
for large momenta differences together with the signifi- 
cant correlation for small momenta differences support 
our assumption. The effect of spin flip in double ioniza- 
tion needs further experimental investigations. 
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